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Abstract 

We describe energy-momentum conservation in relativistic perturbation 
theory in general FRW backgrounds with causal source terms, such as the 
presence of cosmic defect networks. We provide a prescription for a linear 
energy-momentum pseudo-tensor in a curved FRW universe, decomposing 
it using eigenfunctions of the Helmholtz equation. We also construct con- 
served vector densities for the conformal geometry of these spacetimes and 
relate these to our pseudo-tensor, demonstrating the equivalence of these 
two approaches. We also relate these techniques to the role played by 
residual gauge freedom in establishing matching conditions at early phase 
transitions, which we can express in terms of components of our pseudo- 
tensor. This formalism is concise and geometrically sound on both sub- 
and superhorizon scales, thus extending existing work to a physically (and 
numerically) useful context. 
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I. INTRODUCTION 



Considerable challenges are presented by the study of the causal generation of pertur- 
bations seeding large-scale structure formation and anisotropies in the cosmic microwave 
background (CMB) Not only is the analytic treatment of the resulting inhomoge- 
neous evolution equations extremely complicated, but their numerical implementation 
must also circumvent a number of subtle pitfalls before facing up to the severe dynamic 
range limitations of even supercomputer simulations. To date the only quantitative nu- 
merical studies with realistic causal sources, such as cosmic strings [0-^ or other defect 
networks, have been performed in fiat FRW {K = 0) backgrounds Despite positive 

indications about the large-scale structure power spectrum for models with a cosmolog- 
ical constant included [||, these defect networks in fiat cosmologies appear to be unable 
to replicate the observed position of the first acoustic peak in the CMB angular power 
spectrum PJ^-^ — indeed the best results for defects are for 7^ cosmologies |]10 . 



This situation contrasts markedly with the standard inflationary paradigm in which 
reliable predictions about the CMB acoustic peaks are relatively straightforward to make 
and for which there appears to be remarkable accord with recent CMB experiments 



11]. So the question arises as to the relevance and utihty of complicated theoretical 



studies of causal perturbation generation when the simple primordial inflationary models 
appear to suffice. The first motivation is that the confrontation with observation remains 
indecisive, not only because of the significant experimental uncertainties — for example, 
even MAP data will be insufficient to simultaneously constrain both the adiabatic and 



isocurvature infiationary modes, and cosmological parameters [12| — but also because 



good quantitative accuracy has not yet been achieved for the full range of cosmic defect 
theories. For example, even for fiat universes, a subsidiary role for defect networks 
complementing the infiationary power spectrum cannot be excluded. Indeed, claims of 



improved fits in hybrid defect-infiation models |jT3[| are not surprising given the extra 
degrees of freedom available. 

There are a number of mechanisms by which defects can be produced at the end 
of infiation with the appropriate energy scale: Hybrid infiation typically ends through 
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symmetry breaking which generates defects [Q. Phenomenological GUT models have 
been proposed which can produce superheavy strings after inflation W^. 'Preheating' as 
inflation ends is also capable of creating superheavy defects even for low energy inflation 



scales [0]. Given the foundational uncertainties that remain concerning inflation ||T^ and 
the lack of a widely accepted realistic phenomenology, it is only reasonable to continue to 
explore alternative paradigms such as late-time 'causal' generation mechanisms — which 
are not exhausted by defect networks in any case, e.g. 'explosion' |T^, and other source 
models ^ . Moreover, in order to have confidence in cosmological parameter estimation, 
it will be necessary to constrain these alternative models, including the effects of vector 
and tensor modes, and K ^ backgrounds. Here the combination of intrinsic curvature 
and defect sources is particularly interesting. 

Cosmic defects would be expected to contribute to the nonGaussianity of CMB 
anisotropics and the presence or absence of such distinct signatures will provide obser- 
vational tests with which to confront inflation and causal paradigms |]TP[|. A particularly 
exciting prospect is the detection of a CMB polarization signal for which the competing 
models give very different predictions and, indeed, some causal effects can be differenti- 
ated Of course, the discovery of topological defects, which are strongly motivated 
in our high energy physics, would have profound implications for our understanding of 
the early universe. 

Finally, we note that there is now a significant body of work about causal mechanisms 
for structure formation and this has raised a number of interesting issues within general 
relativistic perturbation theory. However, even with most work undertaken in a flat 
FRW background, the number of approaches to the problem almost equals the number 
of papers. A key aim of the present paper, then, is to demonstrate the equivalence of the 
most important of these approaches and to generalise this work to all FRW cosmologies, 
laying the foundations for quantitative studies in curved backgrounds in particular. We 
shall work in the synchronous gauge because of its ubiquity in numerical simulations and 
the greater physical transparency offered by this gauge choice. 

In the literature, treatments of the energy-momentum conservation of individual 
modes in the combined system of gravitational and matter fields have been variously 
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phrased in terms of 'compensation' |]7|pT|, 'integral constraints' p2| , p3| , and the con 



struction of 'pseudo-tensors' to describe the energy and momentum densities and their 
conservation laws |l5|, |2l| , ^ , |2^ , as well as the use of matching conditions across a phase 



transition to set initial conditions The relationship between these notions and the 
initial conditions has been discussed to some extent, in the case of a flat FRW back- 
ground. For general FRW cosmologies, however, the situation is less clear and deeper 
conceptual issues have to be resolved. 

In §0we shall provide a prescription for the construction of a linear energy-momentum 
pseudo-tensor in the K ^ FRW universe. The pseudo-tensor so obtained agrees in 
the fiatspace limit {K — > 0) with the Landau-Lifshitz stress-energy pseudo-tensor r^,^ 
obtained in ref . . We also discuss the philosophy underlying the notion of a pseudo- 
tensor and how its inherently global nature appears to be at odds with theories of local 
causal objects. In § [m] we define energy and momentum with respect to a general FRW 
background manifold. This allows us to calculate conserved vector densities for the con- 
formal geometry of these spacetimes, and to relate them to our pseudo-tensor, giving it 
a local geometrical meaning that is valid on all scales and demonstrating the equivalence 
of the two formalisms. In §|V| we apply the matching condition formalism p2| , pB| to a 
curved universe, and discuss how the residual gauge freedom in the synchronous gauge 
may be exploited to make the pseudo-energy continuous across the phase transition in 
which the defects (or other sources) appear. We also show that we may match the vector 
part of the pseudo-tensor across this transition. We conclude (§0) with a discussion of 
the implications of this work. 

II. A GENERALISED ENERGY-MOMENTUM PSEUDO-TENSOR 

We wish to consider metric perturbations hf^i, about a general FRW spacetime 

ds'^ = a?{'^^u + h^y)dx^dx'' , (1) 

where the comoving background line element in 'conformal-polar' coordinates (r, x, 0, 9) 
is given by 
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(2) 



with the function suik X depending on the spatial curvature K as 



' sinh X , X < , 



= < 



= 0, 



(3) 



X > 0. 



Here, a = a(r) is the scalefactor, for which we can define the conformal Hubble factor 
H — d/a, with dots denoting derivatives with respect to conformal time r. As emphasized 
earlier, we shall adopt the synchronous gauge defined by the choice 



where the trace is given hy h = ha (with the convention throughout that Greek indices 
run from to 3 and Latin from 1 to 3). 

The Einstein equations are given by G^i, = R^y — ^g^uR = i^T^j^^ (with k — SttG), and 
we will separate the energy-momentum tensor T^^ into three parts: 



The background tensor T^i, includes the dark energy of the universe (or cosmological 
constant), while the first order part STf^^y incorporate the stress energy of the radiation 
fluid, baryonic matter, and cold dark matter. The final contribution 0^j, represents the 
stress tensor of an evolving defect network or some other causal sources. This is assumed 
to be small (of order ST^j^^) and 'stiff', that is, its energy and momenta are conserved 
independently of the rest of the matter and radiation in the universe and to lowest order 
its evolution is unaffected by the metric perturbations h^j,,^. 

A. Conceptual discussion and pseudo-tensors in flat {K = 0) FRW spacetimes 

It is interesting also to consider the notion of the energy-momentum tensor of the 
geometry or gravitational field, which we shall denote as i^i^. If it were possible to define 
then we could re-express the perturbed Einstein equations simply as a wave equation for 
h^i, with a source term constructed from the 'complete' energy-momentum tensor, that is. 



= , 



(4) 



(5) 
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the sum T^,y = T^y+t^^. As we shall explain, the linearized Bianchi identities would imply 
that the sum t^j^ is (to linear order) locally conserved r^^ = 0, since it includes all the 
flux densities of matter and gravity (unlike the covariant conservation law T^^ = which 
represents an exchange between matter and gravity). Such motivations for incorporating 
the geometry in a 'complete' energy-momentum tensor t^^j are discussed at considerable 
length in ref. p6| using the example of metric perturbations about Minkowski space. 

Einstein, as well as Landau and Lifshitz, have presented procedures whereby one 
may rewrite the Bianchi identities to obtain quantities that they call energy-momentum 
"pseudo-tensors" . These have some of the above properties, and allow for the calculation 
of various conserved quantities |]2^. Here both t^y and t^^ are quadratic in the connec- 



tion coefficients, so that they are "linear tensors", behaving like tensors under linear 
transformations. 

For a Minkowski space, with 7^^, = 77^,^, a = 1 in (|l]), linearising reveals this procedure 
to be essentially trivial because t^y vanishes to first order. However, for the fiat space 
{K = 0) expanding universe, the time dependence of the scalefactor a in (|l]) introduces 



additional terms at linear order. This has been used by Veerarghavan and Stebbins ||2T 
to define an energy-momentum pseudo-tensor in this case: 

nil 

Too = {STqo + 600) , Toi = 5Tok + 6ofc , 

K 

Tij = 5Tij + Qij (hij - h5ij) . (6) 

Here the components and r*-' defined in @ can be identified as the pseudo- 

energy density U, the pseudo-momentum density S, and the pseudo-stress tensor Vij 
respectively. Using the stress-energy conservation equations (the Bianchi identities) 

rZ = ' (7) 



various suitable choices of evolution variables have then been made: for example, 



This fiat space result can be obtained from a straightforward manipulation of the field 
equations for h^^^ which involves moving any background-dependent terms to the right 
hand side ||^. However, for the generalization to curved spacetime backgrounds we need 
a more rigorous prescription for the energy-momentum pseudo-tensor, as well as the 
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definition of its components in a coordinate system appropriate for practical applications 
— this is the subject of this section. We are also called upon to come to terms with the 
non-local nature of these objects. 

The Landau-Lifshitz construction of r^j^ proceeds by appealing to the principle of 
equivalence, which allows one to choose a normal coordinate system so that the con- 
nection coefficients vanish in the neighbourhood of a point. In a general spacetime, the 
interacting part of the geometry t^^ cannot be made to vanish by this coordinate choice, 
although it then resides only in the second and higher order derivatives of the metric. 
Nevertheless it becomes significant over extended portions of the spacetime and so the 
energy- momentum of the geometry must be understood as global in nature 1^^. This 
fact forbids the existence of a tensor density for the gravitational energy and momenta, 
so that the best that we can actually hope for in terms of local quantities is a "pseudo- 
tensorial" [| object which, suitably integrated over a large region of spacetime, would 
lead to a quantity that is sufficiently gauge invariant for practical purposes. 

However, in causal perturbation theory, we are particularly interested in a distribu- 
tion of small perturbations each of which has associated energy and momentum. These 
objects (such as topological defects, and their associated perturbations) are not well 
modelled, even as a distribution, by quantities that have no meaning except over large 
portions of the spacetime, and one has a rather ad hoc balance between the requirement 
that one consider a sufficiently large volume, and the understanding that effect of the 
distribution of causal objects should average to zero. This is also a major conceptual 
difficulty facing integral constraints for localised perturbations as discussed by Traschen 
et al Fortunately, there exists a formalism |^ in which one can avoid these dif- 

ficulties by defining energy and momentum with respect to a background manifold, so 



These objects are commonly known as "pseudo-tensors" for historical reasons, e.g. Einstein's 
anti-symmetric construction. Here, the nomenclature refers to the fact that they require ad- 
ditional structure — such as a preferred coordinate system/background manifold — on the 
spacetime for their definition |2^, rather than their transformation properties under reflec- 
tions. They are not true tensors, but linear tensors. 
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that one obtains conservation laws and conserved vector densities. We shall apply this 
formalism to the general FRW spacetime in §|T|, thereby providing the results of this 
section with a local geometrical interpretation on all scales. 

B. General FRW {K ^ 0) spacetimes and curvilinear coordinates 

Consider two spacetimes related via a conformal transformation — also known as a 
metric rescaling — of the metric tensor so that 

~9^. = ^9,u , r = ^-'g'" , (8) 

where is a scalar function of the coordinates ^(x'^). A general FRW universe may 
be so rescaled to a stationary (a = 1) FRW universe. Since the non-zero intrinsic 
curvature of a general FRW spacetime manifests itself in the non-vanishing property 
of the background Einstein tensor (even in a stationary spacetime), we shall have to 
separate out the background from the perturbed parts. Moreover, since we wish to 
express perturbations in terms of the Helmholtz decomposition in polar coordinates, we 
shall write all spatial derivatives in terms of the covariant derivative with respect to 
jij, rather than the partial derivatives as previously for the K = case in Cartesian 
coordinates. 

Under (H) the Einstein tensor transforms as 

V = -V'm;!. + l^i^tii^u + ^9i^uVi^<T + Qtiui^^-a > (9) 

where ip^ = (Infi)^^. Now let g^u = CL^ilnu + hf^,^) as in (|l]) with Q = l/a^, so that 
g^ll' = luf + is the metric for observers comoving with the expansion of the universe. 
If we raise the first index, we can make the identification ipQ = — 2?i, ipi = 0. Hence, the 
components of a stress energy "pseudo-tensor" defined by 

r.^ = G'J^ , (10) 

may be written as 
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= -3K + (^a^SG\ + nh) , 

^r'. = -K5'^ + (^a^SG'j-n[h'j-hS'^]) . (11) 

We note that, since metric rescalings (H) preserve its tensorial properties, the r^^ defined 
in ( pup are true tensors in both the stationary and the expanding spacetimes. These 
may be further written as a sum of a background contribution from the unperturbed 
spacetime, and a perturbed part (unhke the K = case for which the background term 
vanishes). Thus, Tj^ = fj^ + Srj^, with the components given by 

Kf^Q = -3K , k6t% = a^6G% + Hh , 

= , /t^rO = a^G^^ , (12) 

fi:f^^. = -K6i , /€5r^ . = a^6&^ - n[h' ^ - hS}] . 

Since the r'^j^ are precisely the Einstein tensor (divided by n) in the conformally related 
stationary spacetime (jfj^i,, they must satisfy the Bianchi identities there. Hence, we know 
that 

Dot'o + DjT\ = , b,T\ + b^T\ = , (13) 

where, denotes covariant differentiation with respect to the stationary 4-metric g^,^. 
Now, using the connections and ([T^P; and working to first order, we may rewrite (|T3|) as 



^^% + Kh- = 0, (14) 



where the bar denotes the covariant derivative with respect to the 3-metric 7jj, and the 
— ^h term is implicit in the covariant derivative DjT\. 

This manner of rewriting the Einstein equations clearly reduces to that of [21| — see 
equation (H) — for i^' = 0, where r^^^ = 5t'^^. However, the equations ( [T^ ) obeyed by the 
5r^j^ are more complicated than (j^) because the non-zero intrinsic curvature manifests 
as a non-vanishing background Einstein tensor, which appears in the Bianchi identities 
for the full spacetime. 
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If we exploit the fact that the vector Pq — Kilhng in g — with components Sq 
may be multiphed with itself to form the (reducible) Kilhng tensor —6q6'^, then we may 
add +Kh/n times this tensor to dr^^^ without disturbing the tensorial properties of the 
perturbed part of the r^^. This amounts to a redefinition of the 00-component only. 
Henceforth we shall consider 5r^^ to be redefined in this fashion so that 

K5r° = /«5r° o-Kh = a^6G^Q + Hh-Kh. (15) 

The new 5r'^^ will then satisfy the concise equations 

^Ao + ^A|. = 0, (16) 
5A,o + Kb- = . (17) 

We may justify this redefinition by noting that the 00-component so obtained is precisely 
the definition of energy (up to a factor a?^) obtained from the conformal Killing vector 
Pq in the following section. Furthermore, since the volume element d^ = dV + dV^ert 
where d^ = a^^dydOdct), divert = {h/2)dV and 7 is the determinant of the spatial 
3-metric, we may interpret the —Kh term as representing the alteration to the fiat space 
energy due to the effect of intrinsic curvature on the volume element. 

Expressing the conservation properties of a general FRW cosmology in this fashion 
is particularly useful, as it produces equations phrased in terms of the spatial covariant 
derivative, which is precisely the language used to express the properties of the Helmholtz 
eigenf unctions Q^'"^ commonly used to describe perturbations in such cosmologies. 

C. The Helmholtz decomposed pseudo tensor 

For perturbations over a curved FRW background, we can no longer make use of 
standard Fourier expansions. Instead, it is usual to employ the Helmholtz decomposition 
using the linearly independent eigenfunctions of the Laplacian in polar coordinates (see 
ref. [0). We expand all perturbation quantities in terms of the eigenfunctions Q*-'"\ 
which are the scalar (m = 0), vector (m = ±1) and tensor (m = ±2) solutions to the 
Helmholtz equation 
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V2qM = yiQ(™) = -fc2QM , (18) 

where the generahsed wavenumber g, and its normahsed equivalent (3 are related to k via 
= + (|m| + 1)K, P = q/ ^J\K\ and the eigentensor has \m\ suppressed indices (equal 
to the rank of the perturbation). The divergenceless and transverse-traceless conditions 
for the vector and tensor modes are expressed via Qi^^'"* = and Y''Qit'^^ = Qif = 0. 



Auxiliary vector and tensor modes may be constructed as follows: 



Ql 



(0) 



(0) 



Q 



(±1) 



)(±i) 



(19) 



The spectra for flat and open universes {K < 0) are continuous and complete for (3 > 0. 
For the K > case, the spectrum is discrete because of the existence of periodic boundary 
conditons. For scalar perturbations, we then have (3 = 3,4,5,... since the (3 = 1,2 
modes are pure gauge [^. Using this decomposition the metric perturbation may be 
decomposed as 



h,, = 2 / d/i(/5) /iL7.,g^°^ + hrQf/ + h'^^Qif + h'f'^Ql;'^ + h'^^^Qf.^ + h'a'^Q 



,(2)^(2) 



f-2)^(-2) 



(20) 



where and hj- represent two 'longitudinal' and 'transverse' scalar degrees of freedom, 
hy^ two vector modes and h^"^ two tensor modes. As well as the transform over the 
'radial' coordinate /3, there is an implicit sum over indices im which label the spherical 
harmonics encoding the angular dependence. 



Decomposing the energy-momentum pseudo-tensor ( |15D in this fashion, we have 
= / d/i(/3)rsg(°) , 



'0 

6t? 



5r] 



J d/.(/3) [2 {ra}Q^^ + rrQ^^ + rff Q^^^/ + ri^'^Q^'f ^ 



T, 



'Q 



(21) 



where the tjv and terms are the 'induced-vector' and 'induced-tensor' modes asso- 



ciated with Qf'^ and Qif^^ auxiliary modes. The quantities ts, Ty^^\ tl, tt and t^^'^^ 
are defined as 
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KTs 



hL + 



KTy 



(±1) 



K_ 

2K 

w 



hi 



-Ka' 



[pf5f + Ps]+mhL-6KhL 



1 



{Pf+Pf) 



KT, 



(±2) 
G 



(±2) 



1^2 



;{±2) 



7 a 



5pf + 
r(0) 



(22) 



(±2) ^ n(±2) 



2nh 



(±2) 
G 



In the second equahty for each of the above equations we have made use of decompositions 
similar to (|2lD for the fluid (subscript /) and source (subscript s) terms, so as to 
write the pseudo-tensor in terms of these variables 0]. The equations (pJ]), ([TtI) yield 
four equations for the remaining four variables: 



TlV 



'IT 



11 
k 



kf, 



IV 



2e 



rL + 2 



K 



Tt 



-2k[P - 2K]-^4 



;(±i) 



(23) 



where we have used the first equation to obtain the final equality in the second. We 
observe that we have six independent quantities: ts and one of tl, tt for the scalars. 



Ty^^^ for the vectors, and Tq^'^' for the tensors. (Note that tl is defined as the spatial 
trace: 6tlQ^°^ = r/.) 

Finally, we comment on the relation of our pseudo-tensor ( p2D to an alternative defi- 
nition given by Uzan et al. in ref. p5|. For perturbations over a curved {K ^ 0) FRW 
universe, there exist several possible (ad hoc) generalisations of the Landau-Lifshitz 
pseudo-tensor, depending upon the particular choice of initial conditions and the man- 
ner in which one removes the residual spatial gauge freedom present in the synchronous 
gauge (see later in §^V]). In ref. p5| matching conditions were used (as an interesting 
aside) to define the tq^ components of the pseudo-tensor as 



-(±2) 



KT, 



UDT 



00 



= V7 



k5Too + fi:6oo + Kh 



Hh 



KT 



UDT 



Ok 



K5n. - 2KdkE 



(24) 



where ^qToo = dkTok, and h, E, h correspond to the formalism of this paper as: h = 
6//i(/3)/iLQ(°), -AE = J ij{(3)hTQ^^\ and h- = h-2AE, A = D^A- 

Apart from providing a prescription for all the components of the energy-momentum 
pseudo-tensor (and in a more elegant decomposition), our definition ( p2D extends and 
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improves upon that proposed in ref. ||2^ on two counts. Firstly, p^) was only given a 



geometrical interpretation on superhorizon scales. The Kh term in the (redefined) 5t^ 



component in ([15| ) replaces a Kh~ term in their definition (|2^ , where their variable 
h~ = h — h"^ is the sum Q{hL + hT/3). The two definitions agree in the superhorizon 
limit, in which case h h~, but our definition (^) and its physical interpretation are 
also valid on subhorizon scales. 

Secondly, there are the limitations inherent in the manner in which the pseudo-tensor 
is defined in ref. Unlike (JT^) the perturbed and background parts of the pseudo- 

tensor are not distinguished. Moreover, their quantity r^i is defined via a conservation 
equation, so that the pure divergenceless part Ty^^\ removed by the derivative in the 
equation (|T6[) is not specified. We shall show (in § [rVD that this part can be recovered 
as a vector quantity to be matched across the transition. Finally, the definition of 
Too = — a^r°Q in (^) and reference differs by a factor = from our r°Q, so 



that it is related (on superhorizon scales only) to the one conserved current Jp^, whereas 
all components of our pseudo-tensor ([T^) can be related to the four conserved currents 



defined in the next section (§p] 



D. Relation to the superhorizon growing modes 

The pseudo-energy br^ (or ts) obtained in this section may be simply related to 
the coefficient of the superhorizon growing modes for the CDM density perturbation 
5c in the radiation- and matter-dominated well as in the curvature-dominated 

epoch. Assuming adiabatic perturbations and ignoring the source terms in a two fiuid 
radiation plus CDM model, it is well known that the CDM density perturbation obeys 
the equations: 



3 r 



4 + H5c - 2 



A = , fie = 1 , fir = 



In both the radiation and matter eras, there exists a superhorizon growing mode pro- 
portional to r^, while in the curvature-dominated regime, this becomes a constant term. 
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If we let the coefficient of this mode be A, then we find that kts ~ ~8A in the radiation 
era, kts ~ —20A in the matter era, and kts ~ 2KA in the curvature regime. Thus, our 
generahsed pseudo-energy essentially tracks the growing mode of the density perturba- 
tion. This is a useful property for numerical simulations (as discussed for example in 
1^), since we can replace 6c with ts, thus avoiding the possibility of spurious growing 
modes sourced by numerical errors. We shall further investigate the inclusion of the 



pseudo-tensor in numerical evolution schemes elsewhere |32 



III. FRW CONFORMAL GEOMETRY AND CONSERVED CURRENTS 



The energy, momentum and their conservation laws for one spacetime may be defined 



with respect to another manifold in an inherently local manner In the context 

of perturbation theory, we already have a background, and it seems logical to employ 
this approach. However, this is not a very compact form of expressing the desired 
conservation laws which, unlike the pseudo-tensor of the previous section, are not phrased 
in terms of the spatial covariant derivative with respect to 7jj, making it incompatible 
with the decomposition of perturbation quantities with respect to eigenfunctions of the 
Laplacian. Here we shall calculate the conserved vector densities for the conformal 
geometry of a general FRW spacetime, and relate these to our pseudo-tensor, giving it a 
geometrical meaning that is valid on all scales and demonstrating the equivalence of the 
two formalisms. 

A. Conserved currents with respect to a FRW background 

The longstanding problem of defining energy, momentum and angular momentum for 



general relativistic perturbations has been considered by Katz et al I^Oj. They provide 
a general formalism by which one can define, for an arbitrary spacetime (M, g^j^) con- 
taining perturbations and any vector ^, conserved vector densities /^(O with respect to 
a background (M, g^^p) and a mapping between M and M. Here, and hereafter, a caret 
shall denote multiplication by = ci^^/l where g = detg^j^i, and 7 = det7^,^. Although 
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one may use any vector ^, it is useful to choose ^ as the conformal Kilhng vectors of the 
background spacetime, so as to exploit its symmetry properties. 

In general, the choice of a particular background is free. However, it makes sense to 
either choose simple backgrounds possessing maximal symmetry or to choose as a back- 
ground one that is already commonly in use in cosmology such as an unperturbed FRW 
spacetime. Conceptually, one might desire a background possessing a maximal Killing 
geometry (spanned by 10 linearly independent Killing vectors), so as to immediately 
generate Noether conserved quantities and currents. Of the general FRW spacetimes, 
only de Sitter spacetime has this property. The implications of the de Sitter Killing 
geometry have been investigated by several authors |^ and it allows for a clear rela- 



tion to Traschen's integral constraints ||23|. We shall demonstrate that the choice of a 
FRW background is not only quite tractable (despite the complications introduced by 
the use of a conformal rather than pure Killing geometry), but also allows for a clear 
relation between the conserved vectors and our 6t^^, which is valid on both sub- and 
superhorizon scales. 

The details of the construction of the conserved vector densities I'^ associated with a 
conformal Killing vector ^ shall be omitted. The general formalism is given in and 



the details of the construction of a relation between the 00-component of a pseudo energy 
momentum tensor and the conserved vector density associated with just the conformal 
Killing vector normal to a constant time hypersurface may be found in [^. We shall 
simply quote those results required for the current analysis: for each conformal Killing 
vector ^ we may define a vector density = yJ—gI^{C) by 

«:F(0 = 5G"^,r + ^^r + /^c^ (25) 

and 

A^i" = \ {R\5% - R%6^,) h^^C = ^[n-n'-K] (26) 

8Ata\^ = {hg^'P - h^P) - Dp {hg^"" - h^'f) (27) 

where we have substituted for the background terms in (|26|), Z = g^^'^Z^.v, and Z 



C^gf_iu = 2ipg^^. Here, £ denotes the Lie derivative, and ip is the conformal factor for ^ 
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with respect to g^,y, so that = for ^ Kilhng. The vector density so constructed will 
satisfy: 







(2^ 



where we have used the result = {y/—gV^)^^/ \f--g for an arbitrary vector V in 



the first equality of (0), and the FRW connection coefficients in the last. Here, and 
elsewhere, we have used the subscript ^ to denote that the conserved vector so labelled 
is generated by the vector ^. 



B. Relations between the 6t^^ and the 

Any conformally flat spacetime will admit a maximal conformal Lie algebra spanned 
by 15 linearly independent conformal Killing vectors. For the general FRW metric, these 
were obtained in ref. in the coordinates (r, x, ?/, z). In principle, each of the 15 vectors 



will generate a conserved vector with four components, and one conservation equation, 
yielding at least 45 components. Given that the symmetric pseudo-tensor has only 10 
linearly independent components, of which 4 are removed by the Bianchi equations (|T6|) 
and ([17|) , there is clearly a considerable redundancy in the information contained in the 
set of all the vector densities obtained using the FRW conformal geometry. Since we 
wish to relate these conserved currents to our pseudo-tensor, our choice of vectors is 
guided by the desire to keep ^ simple (so that 5G^^^^ may be simply related to the 
^'^^lu), and for the vectors to pick out different components t^^. We shall therefore be 
particularly concerned with: the conformal Killing vector Pq normal to constant time 
hypersurfaces with conformal factor ipp^ = Ti; the angular Killing vectors M12 and 
M23; and the generalised isotropic conformal Killing vector H which has the conformal 
factor ipH = cosii: x[7Yn(r) + n'(r)]. In the coordinates {t,x,0,4>), these vectors have 
components: 

Po^ = (1,0,0,0), Mr2 = (0,0,0,l), 

M23 = (0, 0, — sin(/), — cot 6^ cos (/)) , if^ = (cos/^ x "-('*")) siui^ X n'(r), 0, 0) . 

Here, sinx x is defined in (^, while coskX = {coshx, 1, cosx}; and n(r) = {coshr, r, 

cost} for K < 0, i^T = and K > respectively. 
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These conformal vectors reduce to Killing vectors under special conditions on the 
scale factor: for a flat K = FRW spacetime, the vector Pq is Killing if a{t) = C 
where C is some constant so that we have the stationary Einstein spacetime; and H is 
Killing if a(t) = Cexp(— t/C) so that we have a de-Sitter background. In the case of the 
K spacetimes, Pq is Killing if a{t) = C; and H is Killing if a{t) = C/h{T), where 
/i(r) = {cost, coshr} for K = {—1, 1} respectively. 

Using ( p9D in (|25| ) we obtain the following conserved vector densities which relate 
directly to our pseudo-tensor (5r/f given in 



'^-'M23 



sin k5t^2 ~ cot 9 cos kSt^^ 



sin 6 k,6t 9 — cot 6 cos d) k6t o 



cot 6 cos (f)(h''^ - S'^^iij 



sm 



(30) 
(31) 



(32) 



valid for all FRW spacetimes, as well as 



^ K, (5t°q + Kh) cosh T cosh x + kSt^i sinh r sinh x + ^ sinh r cosh x 
K(5r°Qr + (5r°;^r + h 

K (5t°q + Kh) cos r cos x — i^St'^^ sin r sin x ~ h sin r cos x 



K <0 
K = 



K>0 



'^/7 k6t^q cosh X cosh r + K(5r'^-^ sinh x sinh r + 7i sinh x sinh r (^/i'^-^ — 5*^]^/;, 
_ ^fci-j gjj^j^ ^ _j_ ^j^ki _ ^|fc^ cosh x) ('^^ cosh r + sinh r 



H 



k^t^'q cos X cos r — /t(5r\ sin x sin t — Ti sin x sin r (^/i'^;^ — 5^^;^/;, 
+ i — {h'^^^ — h^^) sinx + {h^^ 1/ — /i''^) cosx) (Ticosr — sinr 



is: < 

JsT = 



(33) 



where we have used (|T5|) in (^) for each of Pq, M12, M23, and H. 
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C. Alternative derivation of Jr^^^ from the I^'s 

The constraint equations (0) and (p^ satisfied by the energy-momentum pseudo- 
tensor are encoded in the vector density equations (|28|). For ^ = Pq, the equation (^81) 
yields (||); for ^ = M12 it yields (0) with i = 3; for ^ = M23 it yields (0) for i = 2,3 
in the following linear combination: 

- sin (5r°2,o + ^r'^sifc) - cot 9 cos (5r°3^o + 6t%^,,) = , 
while for ^ = H we obtain (^) and {\l7\ j for i = 1, in the combination 



cosh r cosh x t°o,o + '^\\k + sinh r sinh x 



' 1,0 ' l|fc 







for the K < case, and similarly for 7^ > 0. 



Note that since /^^ = <^==^ /^.^ = 0, the identification of the components of the 
perturbed part of the pseudo-tensor as being proportional to the components I'^ leads 
one to expect a conservation law of the form given in (|T6|) and {\i7\). The presence of 
terms in (pO|^3D other than the 6t^^ accounts for the difference between the general 
covariant derivative on the FRW spacetime, on the one hand, and the spatial covariant 
derivative and temporal partial derivative, on the other. Hence, we see that given the 
for the conformal geometry {^} of the background FRW spacetime, we could construct 



the perturbed pseudo-tensor directly using (|30|-p^) and the final equation of (|28|): the 



two formalisms are equivalent. The results of this section also demonstrate that the 
use of a FRW spacetime as the background manifold has the effect of removing the 
background energy and momentum: there do not appear any contributions from the f^^^ 
in the vector densities I^. 

Approaching the 5r'^^ from this point of view also lends weight to the (apparently) 
ad hoc inclusion of the Kh/n term into the perturbed pseudo-energy 6t^q as defined in 
(pTsP because it is this redefined quantity that appears in the conserved (energy) vector 
density associated with the conformal Killing vector Pq. This is not surprising, as the 
isometry described by the Killing vector Pq in the spacetime (M, g) is not entirely lost 
as we go to the spacetime (M, g), where Pq is a conformal Killing vector. It is preserved 
in the evolution space TZ x TM — where TZ accounts for the affine parametrization of the 
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geodesies, and TM is the tangent bundle — by the appearance of an irreducible Killing 
tensor if = —0^5^51 + a^5^, related to the reducible Killing tensor L^^ = —5^5^ + 5^ 



in the (M, g) spacetime |^ . As this last tensor is reducible (a sum of products of 
the Killing vector Pq and the metric, with constant coefficients), it encodes the same 
information as the Killing vector itself. Thus, we may expect there to be an "energy 
isometry" associated with the tensor 5^5^, which we used in § |1IB| . 

The vector densities of this section provide a consistent definition of energy and mo- 
mentum with respect to a FRW background and, as we have just shown, the identification 
of the quantities 5t^^^ (including the curvature term in the 00-component) leads naturally 
to a concise and algebraically useful conservation law, phrased as a differential equation. 

IV. ENERGY-MOMENTUM PSEUDO-TENSORS AND MATCHING 

CONDITIONS 

We wish to consider the emergence of a topological defect network (or other causal 
sources) at some stage in cosmic history, that is, the time when defects 'switch on' 
and are carved out of the background energy density during a phase transition. This 
process sets the initial conditions for all the perturbation variables prior to their sourced 
evolution, a state we must specify if we are to perform realistic numerical simulations. 
It is common to assume that any phase transition at which defects will appear will 
take less than one Hubble time, so it will be effectively 'instantaneous' for all modes 
larger than the horizon at the time of the transition. Matching conditions have then 
been found to relate the resulting perturbation variables on superhorizon scales to their 



prior unperturbed state in a 'sourceless' universe While this approach will apply in 
many physical situations, there are circumstances in which it may not, such as hybrid 
scenarios with mixed perturbation mechanisms or late-time phase transitions in which 
subhorizon modes might be important. Here, we have already defined a generalized 
energy-momentum pseudo-tensor applying to both sub- and superhorizon scales which 
should prove useful for this wider class of scenarios. We shall now demonstrate, in an 
appropriate synchronous gauge, that its components can be used to specify the matching 
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conditions valid for all lengthscales in a defect-forming transition. 



A. Matching conditions on a constant energy density surface 

If the phase transition appears instantaneous for a given mode, we need only to match 
the geometric and matter variables on the spacelike hypersurface surface S, described 
by the equation 

p(x^) = pQ + 5p = const. , (34) 



where, up to a small perturbation, we have assumed homogeneity on either side of S p2 
Prior to the phase transition, the perfectly homogeneous and isotropic 'perturbation' may 
always be absorbed into a redefinition of the (continuous) scale factor. In a simple model 
without surface layers {i.e. ignoring the internal structure of the phase transition), 
the standard procedure used to match the geometric and matter variables is to insist 
that the induced 3-metric _L^,y and the extrinsic curvature K^^j must be continuous over 
S. This task is simplified if, on either side of the phase transition, one uses the residual 
gauge freedom in the time coordinate r — > f = r + T, with T a non-trivial first order 
scalar function of the coordinates, to transform to a coordinate system in which E is 
defined by the equation f = const, (p = const.), and 6p = 6p + poT = 0. Using the 
Friedman equations the appropriate transformation is therefore specified by 

^ _ Sp _ na^ [p6 + p'] 

po m{'H^ + K){l + u)' ^ ' 

which may be interpreted (at each point in 3-space) as moving the time-slicing for- 
ward/backward so that the surface E is a constant time hypersurface. Here p = up 
is the equation of state for the total fluid, but for the purposes of this paper, we may 
assume that we are in the radiation dominated epoch. 

In setting up this gauge, no use is made of the residual scalar freedom, — > x^ = 
x^ _l_ JD^i^^ in the spatial coordinates. Here D^L = d^L because L is another first order 
scalar function of the coordinates. Note that this new gauge cannot be comoving, as 
this would require that T = 0, and we need this freedom to force the constant time and 
constant energy density surfaces to coincide. 
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B. Matching the scalar modes 



In the gauge described above — denoted by a tilde — the metric is given by: (7^1, = 
Ifiu + hfj^i, where we shall rewrite the spatial metric perturbation as 

hij = 2\iLlij + 2 [OiDj - , (36) 

where hj;^(r, x'^) = / d/i(/5)/ii(r, /9)Q^°^(r, x*^, /3) and similarly for hj'(r, x^). These spa- 
tially dependent variables are used as the physical interpretation of the transformation 
is more transparent, and they facilitate comparisons to existing work [^,^. We shall 
obtain results for the /? dependent quantities later. 
The gauge transformed quantities are given by: 

hoo = hoo + -2 fr + HT) , hoi = hoi + Lu - T\ 



00 = "-00 -I — ^ yi -\- ni J , iLQi = lioi -t- — i\i ■, 

hi = hi + + , liT = hr + L . (37) 

Preservation of synchronicity (/iqo = = h^i) thus provides the form of T and L: 

T=l^, L = g{x') + f{x') I - . (38) 

a J a 

where /, g are functions of the spatial coordinates only. As noted previously, / is com- 
pletely determined by the process of establishing a time-slicing that also has constant 
energy density (at the phase transition). However, g is completely free, and may be 



chosen in such a manner as to simplify equations [0. We shall demonstrate that this 
freedom may be more profitably used to specify gauges (for both K = and 7^ 0) in 
which the energy- momentum pseudo-tensor of §|I| must be continuous across the phase 
transition. 

The vector orthonormal to the constant time hypersurface is given by ra^ = 
(—a, 0, 0, 0), so that the perturbed parts of the induced metric _L^j, and extrinsic curva- 
ture Rf^i, are 

S ±ij= a%j , 6 ±^0= , ^^^^ 
SK^, = = SK\ , 5K^^ = , 
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where we use a(r) ^ + Ti-T], obtained by Taylor expanding about r. Assuming 

that the background is continuous across the phase transition, we need only match the 







6K' 



perturbed parts; i.e. we insist that 
lim,_^o+ [F{tpt + e) - F{tpt - e)]. 

Substituting (|5BD, transforming back to the original gauge and using 

nf , 1 



where [F]± denotes the limit 



dr 

a 



f ( -m' + K){l + u) ^ ^ \ ^ lAf 



DiDj - -7,, A hT + 



3 a 
dr^ 
a ^ 

f 



± 



a 










we find that 
(40) 

(41) 

(42) 

(43) 



± 



Taking the linear combination — 67ix (|4l| ) + GKx (^OD we have 



Ts + 2KUg + Af J 



dr 
a 



2n- 







(44) 



where we have used ( |5B| ) and 

Decomposing with respect to the Helmholtz equation, and noting that the eigenfunc- 
tions separate and are time independent, we obtain 



/..(/?) + ^-y^(/5)-y/(/3) 



dr 
a 







± 



2/ 
a 







J ± 



hT{(3)+g{/3) + fi/3) / — 

a 



dr 
a 

m 



J ± 








(45) 
(46) 
(47) 
(48) 



where we have replaced (^) by (PI). 

There exists an entire class of objects related by gauge transformations to the "pseudo- 
energy" 5r°Q corresponding to different choices for g{x^) in (|4^). Uzan et al make 
use of this freedom to specify 



g = -Ht - f J 



dr 

a 
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which ehminates the matching condition (|^) and yields [T(5^^^]± = 0, refer to (0). 
On superhorizon scales this reduces to a matching on our pseudo-energy: [ts]± = 0. 
However, one is not using the gauge freedom to relate the matching condition to well- 
defined geometrical objects. It would be both more aesthetically appealing and more 
useful if one could employ this freedom to make 6t^q continuous across the transition. 
This is a subtle issue that shall be more fully explored elsewhere |3^, where we discuss 
initial conditions and their consistency with causality. For now, we merely note that, for 
practical purposes in which we wish to describe the onset of defect induced perturbations 
carved out of the background (or inflationary) fluid, compensation between the fluid and 
the source densities implies that we can usually take / to be continuous across the 
transition. In the absence of primordial density perturbations, it will moreover initially 
vanish — see equation (^5]). In this physical context, we may then completely specify the 
gauge by choosing = so that we obtain: 

[hLm± = , [Tsm± = , 

N(/3)]±=0, [hT{f3)]^ = 0. (49) 



C. Matching the vector and tensor modes 

The residual gauge freedom in the vector modes may be expressed as invariance under 
the infinitesimal coordinate transformation ' + L\ where L(r, x'') is a diver- 

genceless 3-vector: DiU = 0. Writing hij = 2 (^hy^ ^-^^^ + hy^ ^-^^^^ for the spatial metric 
perturbation, where liyp {t,x'^) = J dfi{j3)hy^^\T, i3)Q['^^\i3 , x^) , the gauge transformed 
vector quantities are 

h.oi = hoi + Li, h.vi = hyi + Li. (50) 

Preservation of synchronicity implies that Lj = everywhere, so that L is a function of 
the spatial coordinates only. Proceeding as for the scalar perturbations we match the 
induced metric and extrinsic curvature. 

After transforming back to the original gauge, and exploiting the fact that Lj = 
everywhere so that D^^Lj-) = is certainly true on the hypersurface, we obtain 
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(51) 
(52) 



Helmholtz decomposing and assuming that /3 modes separate, (^2D is equivalent to 



as the m = ±1 contributions are hnearly independent. Hence, we find that 



(53) 







(54) 



This is precisely the divergenceless part of which is not obtainable by integrating 
the conservation equation (pSf), unlike the induced vector mode r]y (constructed from 
scalars). Using (^2|) we see that, the matching condition ( p^ implies a 'compensation' 
between the source and fiuid vorticities. We shall investigate this phenomenon further 
in the context of establishing consistent initial conditions in ref. [^. The remaining 
equation (^) may be written as [hy^^\l3)]± = by means of an appropriate specification 
{Li = 0) of the residual gauge freedom in the vector mode. 

For the gauge invariant tensors, we find that the (Helmholtz decomposed) tensor 
metric quantities are constrained to be continuous across the transition: 



k 



(±2) 
G 







,(±2) 



± 



(55) 



There is, however, no residual freedom in the tensor modes, so these matching conditions 
cannot completely constrain the continuity properties of the pure tensor contribution 



Tq . It may also be permissible to insist that Tq^^' = initially, although this is not 
mandated by our results. 



(±2) 



V. DISCUSSION 



In this paper we have considered definitions and conservation laws for quantities that 
may be used to define energy, momentum and the stresses, which are of relevance to 
setting the initial conditions for and/or constraining the evolution of numerical simula- 
tions. To this end, we have constructed an energy-momentum pseudo-tensor for FRW 
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cosmologies with non-zero curvature and have generated conserved vector densities using 
the conformal geometry of a general FRW background manifold. We showed that these 
two formalisms are equivalent so that the pseudo-tensor components are geometrically 
well-defined objects on all scales. These results hold in the presence of a non-zero cosmo- 
logical constant, as all the quantities discussed here are purely geometrical constructs, 
describing the symmetry properties of the FRW spacetime. This pseudo-tensor is likely 
to be a useful tool for detailed investigations of causal models in curved FRW universes, 
particularly for hybrid models with mixed primordial and causal perturbations. We have 
phrased these results in terms of the commonly employed Helmholtz decomposition with 
respect to the eigenfunctions of the Laplacian. 

We considered an instantaneous phase transition early in the universe as a first ap- 
proximation to a model for the defects "switching on" , and employed a gauge in which 
constant energy and constant time surfaces coincide. Matching conditions then imply 
that there exists an entire class of objects which are continuous across the transition 
and are related by gauge transformation to our pseudo-tensor components. The no- 
tion of compensation together with a particular gauge specification removes this redun- 
dancy such that the ts (pseudo-energy) and Ty^^^ ( diver genceless vector) components 
of our generalised pseudo-tensor have this property. For a universe which was unper- 
turbed (and hence homogeneous and isotropic) prior to the transition, we may then take 
Ts = = Ty^^^ as natural initial conditions. This result is true on all scales. In a 



subsequent paper we shall establish with more rigour the effect of causality on the 
superhorizon behaviour of the energy and momentum in general FRW cosmologies, as 
well as the implications for setting the initial conditions. 
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